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We investigate the effect of a ferromagnetic dot on a thin- 
film superconductor. We use a real-space method to solve 
the linearized Ginzburg-Landau equation in order to find the 
upper critical field, H C3 . We show that H C3 is crucially de- 
pendent on dot composition and geometry, and may be sig- 
nificantly greater than H C2 . H C3 is maximally enhanced when 
(1) the dot saturation magnetization is large, (2) the ratio of 
dot thickness to dot diameter is of order one, and (3) the dot 
thickness is large. 
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I. INTRODUCTION 

Recent experiments J^-Q] involving thin-film supercon- 
ductors have investigated the effects of nanosizc artificial 
pinning centers in the form of ferromagnetic dots. Mag- 
netic dots with diameters on the order of 200 nm and 
thickness on the order of 40 nm were fabricated on a su- 
perconducting Nb film by electron beam lithography. It 
was found that a regular array of magnetic dots can dra- 
matically influence the transport properties in the pres- 
ence of an applied magnetic field. In particular, the re- 
sistivity displayed minima at "matching fields" in which 
the number of flux quanta per dot was an integer. Such 
effects arc known to occur in regular arrays of empty 
holes in thin-film superconductors Q. These effects are 
believed to arise because empty holes can form effective 
pinning centers for multiple flux quanta vortices Q , lead- 
ing to particularly stable configurations at the matching 
fields. The strong pinning leading to such multiple vor- 
tices arises due to an enhancement of the order parame- 
ter near the edges of the hole, in a manner analogous to 
surface superconductivity §]-§]. 

However, the physical situation for magnetic dots is 
significantly different. For empty holes, the supercon- 
ducting order parameter must have a vanishing derivative 
at the vacuum/ superconductor interface. In a mag- 
netic field, the states satisfying the linearized Ginzburg- 
Landau equation and this boundary condition turn out 
to have a maximum just inside the superconductor, lead- 
ing to a magnetic field H at which the superconducting 
order parameter may be non-zero near the surface while 
it vanishes in the bulk of the sample (i.e.,H > H C2 ). This 
is very similar to surface superconductivity and the 
maximum field at which superconducting order survives 
in the empty hole system is a direct analog of H C3 MM ■ 



For magnetic dots, the strong field present inside the 
ferromagnet supresses the superconducting order param- 
eter, and in such situations it is appropriate to adopt 
a boundary condition in which the order parameter it- 
self vanishes. This spoils the effect that leads to sur- 
face superconductivity, and it is not at first obvious why 
magnetic dots should support the relatively large super- 
currents associated with multiple vortices. However, the 
problem of a ferromagnetic dot in a superconducting thin 
film has a dimension not present in the empty hole ana- 
log: the magnetization and fringing magnetic field of the 
dot itself. (Throughout this work, we will consider only 
dots small enough to be treated as single domain ferro- 
magnets.) For dots with diameter 2R larger than their 
height t (see Fig. 1) ||, shape anisotropy dictates that 
the magnetization of the dot will lie in the plane of the 
superconductor Jlcfl in the absence of any external field. 
The application of a perpendicular magnetic field tilts the 
magnetization out of the plane of the dot, introducing a 
component to the dot fringing field that partially cancels 
the external applied flux passing through the supercon- 
ducting film. This introduces a region just outside the 
dot in which the net field intensity is smaller than the ap- 
plied field, allowing an enhancement of superconducting 
order. 
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FIG. 1. Magnetic dot with a radius R and a thickness 
t(< 2R) at the center of a thin-film superconductor. An ex- 
ternal field Hqz is applied through the sample, where z is the 
direction parallel to the normal of the sample plane. 

To demonstrate this effect, in this work we study the 
analog of H Ca in the presence of a single ferromagnetic 
dot with a diameter greater than its height (2R > t). To 
do this, we solve the linearized Ginzburg-Landau (GL) 
equation using a real-space method to be described be- 
low. The resulting equations specify a maximum mag- 
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netic field H C3 (l) at which a non-zero superconducting 
order parameter may be present for each value of the 
vorticity I. A typical example of our results is shown 
in Fig. 2. The form of this figure is easily understood 
if one keeps in mind the analogy between the linearized 
GL equation and the problem of electrons in a magnetic 
field In this analogy the vorticity I plays the role 

of angular momentum, and it is well known that, in the 
absence of a dot, the lowest-lying single particle states 
of a given I in a magnetic field are localized near a ra- 
dius Ri = a m V2l, where a m = yj is the magnetic 
length, e* the charge of the carriers, and Hq is the applied 
field. While the presence of the dot and the appropriate 
boundary condition changes the precise relation between 
I and i?;, it nevertheless remains generally true that Ri 
increases with I, Thus, for small values of I, H C3 (l) is 
supressed due to the boundary condition on the order 
parameter, whereas for large I, H C3 (l) — > H C2 , the value 
one expects in the absence of the dot. For intermediate 
values of I, one generically sees a peak in H C3 (l), due to 
the fringing field effect described above. 
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FIG. 2. Magnetic field versus vorticity curve for a Nb su- 
perconducting film with a Ni dot (i?=100 nm and £=40nm) at 
T=8.2 K. For large I, the field is unaffected by the presence 
of the dot, giving H C2 = 915.012 G. The maximal H C3 (= 
929.379 G) occurs at 1=8. The inset shows that the difference 
between H c and H C2 is maximal at 2=8. 

This peak value of H C3 (l) gives the maximum applied 
field in which the thin film may sustain superconducting 
order, and is the analog of H C3 . That it occurs at a finite 
value of / indicates that it is indeed true that magnetic 
dots support and presumably pin multiple vortices. It is 
interesting to note that for appropriate parameters, this 
peak may become quite pronounced, and that it may 
exceed the value of H C3 = 1.695ff C2 that occurs for a 
simple infinite surface and represents the maximum value 
possible in an empty hole ||. Such large values of H C3 
occur when (1) the saturation magnetization of the dot 



is large, (2) the height of the magnetic dot t is large, and 
(3) the dot diameter 2R is close to the height t. 

This article is organized as follows. In Section II we 
describe our model of a ferromagnetic dot in a thin su- 
perconducting film in detail and present the correspond- 
ing GL equation for the system. Section III outlines our 
method for solving the equation, and in Section IV we 
present our results. We conclude in Section V with a 
summary. 

II. MODEL OF A FERROMAGNETIC DOT IN A 
THIN-FILM SUPERCONDUCTOR 

Consider a thin-film superconductor with a small mag- 
netic dot at its center (Fig. 1). The radius of the dot is 
R, and the thickness is t. We assume that the dot is 
small enough so that its magnetization density is uni- 
form throughout the dot; i.e., there are no domains. We 
wish to compute the largest magnetic field H C3 for which 
the order parameter is non-vanishing inside the super- 
conductor. 

The magnetization of the ferromagnet is maximal, but 
its direction may vary. The orientational energy in any 
given direction can be conveniently estimated if we ap- 
proximate its cylindrical shape as an ellipsoid whose 
scmimajor axis length is c(= R) and semiminor axis 
length is a(= t/2), as shown in Fig. 3. 
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FIG. 3. (a) Oblate spheroid with semimajor axis c and 
semiminor axis a. (b) z//a. The angle between M s and the 
sample plane is 0. 

The energy for a given orientation of M (the dipole 
moment of the ferromagnet) in an ellipsoid is [ ficfl 
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i*M = ^[(M s cos9) 2 N c + (M s sin0) 2 A Q ], (1) 

where f2 is the volume of ellipsoid, M s is the saturation 
magnetization of the ferromagnet, and N a and N c are the 
demagnetizing factors along a and c. 9 is the angle be- 
tween M and the sample plane. Since N a > N c , sin 9 = 
in the absence of an applied magnetic field and the mag- 
netization is perpendicular to the a-axis, i.e., anywhere 
in the sample plane. 

When there is an external field Hqz (i//a), we must 
add — VLM ■ H z to the dot energy. Therefore, the total 
magnetization energy of the dot is 



-QM S H sin 6 + F^(9). 



(2) 



For a thin-film superconductor, the free energy may be 
written as 

F S =F N + [ d 2 r[a\^\ 2 + i/3|^| 4 + 

Jr>R 1 
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Since the superconductor is thin, we ignore any fields 
produced by supercurrents in the film. The vector po- 
tential A then is the sum of the vector potential due to 
the magnetic dot (A m ) and the vector potential due to 
the external field (Aq): A = A m + Aq. In a uniform ex- 
ternal field Hqz, Aq, in cylindrical coordinates, is given 
by 



A = ]^H r(p. 



(4) 



Furthermore, for a thin-film superconductor, magnetic 
field components parallel to the film have no effect on 
the order parameter [pi. Thus, we may write 



A M (r) = M s tsm9a (—)ip, 
H 



(5) 
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—QM S H + M 2 n(N a - N c ) sin (9. 



(8) 



To find H C3 , we consider \ip\ << 1 and so drop terms 
involving ip. Then we have 
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M s (N a -N c ) 



(9) 



After minimizing the free energy with respect to tp, we 
obtain 



(10) 



where a oc T — T c . Eq. (10) is the linearized Ginzburg- 
Landau equation. It should be kept in mind that al- 
though the equation is linear in tp, it is non-linear in sin 9. 
This means our solution of Eq. (10) requires a level of 
self-consistency usually absent in solving the linearized 
GL equation, which we describe below. 

Eq. (7) essentially defines our model of the film-dot 
system, and Eqs. (9) and (10) are what need to be solved 
to obtain the critical field of the system. In the next 
section we present some details describing how this is 
done. 



III. NUMERICAL SOLUTION OF LINEARIZED 
GINZBURG-LANDAU EQUATION 

We need to find the largest value of Ho (i.e., H Cs ) for 
which the eigenvalue equation (Eq. 10) has a non-trivial 
solution with the boundary condition 



ip(r = R) = 0. 



(11) 



where 
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which is an elliptical integral [Q. Therefore, the total 
free energy of the system is F = Fs + Fm , or 

F = F N + f d 2 r[a\ij\ 2 + \l3\4>\ i + 

Jr>R * 



1 
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"A 



QM S H sin 9 



^l[(M s cos9) 2 N c + {M s sm9) 2 N a ]. (7) 

We need to minimize the total free energy with respect 
to tp and 9 (or sin0). After minimizing with respect to 
sin 9, we obtain 



Due to circular symmetry, we can write the solution to 
Eq. (10) in the form 



m = /(r)e^, 



(12) 



where I is an integer. The orbital number, Z, corresponds 
to the vorticity of the solution for the superconducting 
order parameter. Next, we scale out length by defining 



P 



(13) 



where the magnetic length a m is 
am = ( J*_)V". 



(14) 



Applying Eq. (12) and Eq. (13), Eq. (10) becomes 
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{--TPT + (" - = ef(p), (15) 

p dp dp p 



where 



and 



7/ % 1 , r sin^ ,a m p . 
A(p) = -p + M s t— r a (^) 



2m\a\ , he 



(16) 



(17) 



The smallest eigenvalue of Eq. ( |l5j ) corresponds to the 
largest ifo for which there is a non-vanishing order pa- 
rameter. Once we find Sq, we can directly write 



H, 



2mc\a\ 
he*Eo 



(18) 



To solve this eigenvalue problem, we used a real space 
method as follows: 

(1) Guess Hq = H C3 , Notice that Ho enters explicitly 
in the vector potential and sin#, and thus cannot 
be scaled out as would be the case for an empty 
hole. 

(2) Define a set of N points p\,p2, /O3, . . . , pn and po = 

R 
a 

(3) Turn the derivatives into differences, thereby trans- 
forming the differential equation into a difference 
equation. 

(4) Set up the column vector 
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and define AX n+1/2 = p n+1 - p n , AX n _ 1/2 = p n - 
Pn-i, turning the differential equation into a matrix 
equation: 
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where 
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(AX„ +1/2 + Al n _ 1/2 ) 3 
1 



/2 



pn(AX n+1 / 2 + AX n _i/ 2 ) ' 



(20) 



Pn = 



(AX n+1/2 + M„_ 1/2 ) 2 

Pn 



(21) 



and 



In - -JX^ T A Y \3^ X n+l/2 + 



Pn(AX n+l/2 + 



(22) 



(5) Diagonalize the matrix to find the lowest eigen- 
value, £oi which gives the highest magnetic field 

jj* 2mc\a\ 

he*€Q 

(6) Finally, one must check if the solution is self- 
consistent. If H* [eo] equals the guessed Ho in step 
(1), then Hq = H C3 . However, if ff*[eo] is not 
equal to the initial guess, then set H* — > H and 
return to (1). The element of self-consistency arises 
because one must determine the orientation of the 
dot magnetization, 9, in the field H C3 . Note that if 
the field H C3 is large enough, then sin 9 = 1 and the 
magnetization is fully parallel to the applied field; 
in this situation Eq. (9) is not appropriate (except 
precisely at H /(M s (N a — N c )) = 1), and we do 
not need to iterate the equations. 



IV. RESULTS 

One of the interesting results that was found in this 
study is that H C3 may be very large for this system. In 
order to examine when this happens, we studied how the 
variables, M s , t and a/c (or N a — N c ) affect H C3 . First, 
we fixed t and a/c to compute H C3 for different values 
of M 8 . Clearly, we expect a large M s to create large 
fringing field, capable of cancelling large applied fields 
near the dot. Presumably this will lead to large values 
of H C3 . Fig. 4 illustrates this for a 160 nm thick dot 
with a/c(= t/(2R)) = 0.9 (or N a = 4tt x 0.361472 and 
N c = 4tt x 0.305689) in a Nb film at T=8.2 K, with 
Ni ferromagnetic dots (M s =509 G). (These parameters 
correspond to those of the magnetic dots studied in Ref. 
EL) For comparison, we also show H C3 for a Dy dot 
(Af s =2920 G) with otherwise the same parameters of the 
system. The enhancement of H C3 for the larger value of 
M s is quite apparent. 
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FIG. 4. Magnetic field versus vorticity curves for a Nb su- 
perconducting film with a dot at T=8.2 K. The dot size is 
fixed: t =160 nm and t/(2R) = 0.9. If the dot material is 
Dy (Af s =2920 G), H C3 = 2080.08 G at Z=33. However, if the 
material is Ni (Af s =509 G), H C3 = 1285.03 G at 1=11. 

We also note in this figure and several that follow that 
there is an apparent jump in H C3 (I) vs. I. If I were a con- 
tinuous variable this would not be a discontinuous jump 
but rather a continuous (albeit sharp) drop in H C3 (l). 
Nevertheless, the sharp behavior is a direct result of the 
nonlinearity of the equations in sin 6. The behavior rep- 
resents a sharp crossover as a function of I in which the 
superconducting order parameters are localized relatively 
close to the dot and ones in which they are further away; 
in the latter case the dot potential is a relatively weak 
perturbation on the result in the absence of a dot. 

Next, we fixed M s and a/c to study the influence of 
dot thickness t on H C3 . For large values of t one again 
expects the field generated by the dot can cancel a rel- 
atively large external field, leading to an enhancement 
of H C3 . One can imagine a situation in which both the 
film and the dot thickness t are the same and are varied 
together. The results reported here will apply provided 
the superconductor is effectively two-dimensional - i.e., 
the coherence length must be larger than t. For many 
interesting materials that may not be possible; however, 
it is quite possible and often appropriate to consider sys- 
tems in which the dot thickness is different than that of 
the film. This is the situation depicted in Fig. 1. Fig. 
5 illustrates H C3 for a Ni dot with a/c=0.9 in a Nb film 
at T=8.2K. As expected, thicker dots (bigger t) indeed 
yield higher values of H C3 . 



FIG. 5. Magnetic field versus vorticity curves for a Nb su- 
perconducting film with a Ni dot at T=8.2 K. The dot thick- 
ness varies from 160 nm to 260 nm, while keeping the value 
of t/(2R) fixed. H C3 increases with t. 

Finally, we fixed M 8 and t and considered the effect 
of a/c on H C3 . The aspect ratio of the dot is relevant 
because it enters into the demagnetizing factors N a and 
N c . Physically, if a/c is close to one, the dot should ap- 
proach a limit in which it is relatively easy to tip the 
magnetization out of the plane p3[ . This again maxi- 
mizes the fringing field available to cancel the applied 
field. Fig. 6 illustrates H C3 for a/c— 0.8 and a/c— 0.9 
(or N a - AT C =1.4860487 and N a - iV c =0.7009889, respec- 
tively), and demonstrates that greater values of a/c (or 
smaller values of N a — N c ) give higher values of H C3 . 
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FIG. 6. Magnetic field versus vorticity curves for a Nb su- 
perconducting film with a Dy dot at T=8.2 K. The dot's 
thickness is fixed at i = 160 nm, but the radius varies: 
t/{2R) = 0.9 (or N a - 7V C =0.7009899) and t/(2R) = 0.8 (or 
N a - JV C =1.4860487). The curve for t/(2R) = 0.9 shows a 
greater H C3 than that for t/(2R) = 0.8. 
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V. SUMMARY 

In this work, we studied H C3 for magnetic dots in a 
thin-film superconductor. To find H C3 , we used a real- 
space method to solve the linearized Ginzburg-Landau 
equation. We showed that the enhancement of the order 
parameter crucially involves the shape anisotropy of the 
magnetic dot and is qualitatively different from that for 
empty holes. We found the enhancement to be maximal 
when (1) the dot saturation magnetization M s is large, 
(2) the dot thickness is large, and (3) the value of t/(2R) 
of the dot is close to f , for which the demagnetizing pa- 
rameters iV a and N c are of comparable magnitude. 
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[13] Note that if a/c is greater than one, then the equilibrium 
direction of the magnetization will be along the cylinder 
axis rather than in the plane. In this case fringing fields 
will be present and disrupt the superconducting order 
even in the absence of an applied field. This situation 
was considered in Ref. Q above, and we do not address 
it in our present work. 
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